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Abstract Whereas penumbral models during the last 15 years have been successful 
in explaining Evershed flows and magnetic field inclination variations in terms of flux 
tubes, the lack of contact between these models and a convective process needed to 
explain the penumbral radiative heat flux has been disturbing. We report on recent 
observational and theoretical evidence that challenge flux tube interpretations and 
conclude that the origin of penumbral filamentary structure is overturning convection. 
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1 Introduction 

Sunspot magnetic fields and dynamics have been studied scientifically for 100 years. 
Despite considerable progress during the last decade, a theoretical framework that 
explains sunspot fine structure, dynamics, magnetic fields and energy balance in a 
consistent manner is only now beginning to emerge. This situation can partly be at- 
tributed to the small horizontal scales associated with sunspot fine structure and the 
relatively poor spatial resolution achieved with spectropolarimetric observations. In ad- 
dition, realistic numerical 3D MHD simulations of sunspots have only recently become 
possible. 

During the last few years, there has been a remarkable improvement in the qual- 
ity and diversity of observational data relevant to the understanding of sunspot fine 
structure, dynamics, magnetic fields and energy balance. In particular, high-spatial 
resolution observations from the Swedish 1-m Solar Telescope (SST) and the Solar 
Optical Telescope (SOT) on Hinode reveal new sunspot structure and flow patterns 
at odds with prevailing interpretations in terms of flux tube models. In addition, the- 
oretical arguments as well as recent 3D MHD simulations of sunspot fine structure 
underline problems of these interpretations and lead to the conclusion that the origin 
of penumbral fine structure is overturning convection. 
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In the present review, we describe recent progress in our understanding of penum- 
bral fine structure and put that in context with existing models. Rather than attempt- 
ing to summarize the extensive literature on penumbrae, we discuss selected key pa- 
pers and attempt to describe their interconnections and to critically review conclusions 
drawn. We also point out connections between observed penumbral fine structure and 
magnetic flux concentrations outside sunspots, such as faculae. We hope to convince 
the reader that the new picture emerging is one of improved consistency as regards 
observations and theory of sunspot penumbrae in particular, but also with respect to 
umbral dots, light bridges and faculae. 



2 Overview of Established Models 

For the past 15 years, the predominant paradigm of penumbral filaments has been based 
on nearly horizontal flux tubes, portrayed and modeled as radially aligned cylinders, 
embedded in a more vertical magnetic field. T hese flux tube interpretations have their 
roots in the work of lMever fc Schmidtl (|l968l ) who proposed to explain Evershed flows 
as siphon flows, originating from a difference in magnetic field strength between the 
two foot points of a flux tube. In the 70's, flux tubes and clusters of flux tubes were 
also established in models of magnetic flux concentrations, surrounded by field-free gas, 
with scales ranging from less than 100 km to that of a large sunspot. 



2.1 Embedded flux tubes 

Of particular import ance in the current s u nspot literature is the uncombed penumbra 



model, proposed by ISolanki fc Montavonl 1 1993h to explain the strongly asymmetric 



Stokes V profiles observed on the limb side penumbra for sunspots away from disk 
center. This model addressed an apparent problem of very strong line-of-sight (LOS) 
gradien ts in the inclination angle of the p enumbral magnetic field, inferred from Stokes 



gradien ts m the inclination angle o t the p enumbral magnetic held, interred trom btokes 
data bv lSanchez Almeida fc Litej ( 19921 ). The large LOS inclination gradients derived 



from this data were (incorrectly, see Sect. 4) interpreted to imply volume currents and 
associate d curvature forces strong enough to completely disru pt static force balance in 
the spot (San chez Almeida fc Lit es 1992; S olanki et al.|[l993) . The uncombed penum- 



bra model avoids this problem, at least partly, by postulating the existence of discrete 
flux tubes, within which the magnetic field is assumed to be homogeneous and there- 
fore current-free. The uncombed model thus 'replaces' smooth inclination gradients 
and (assumed) large volume currents with discontinuous changes at the boundary of 
the flux tube and an associated current sheet. Solanki and Montavon demonstrated 
that a nearly horizontal flux tube, with a strong flow parallel to its magnetic field, can 
explain the observed net circular polarization (NCP) resulting from Stokes V asymme- 
tries of this configuration. Moreover, if such a flux tube is located entirely above the 
photosphere, both its upper and lower boundaries contribute with the same sign to the 
asymmetry of its Stokes V profile, thus enhancing the NCP. 

Various implementatio ns of flux tube mode ls with polarized radiative transfe r 
were later developed by e.g. [ Martinez Pilletl ll 200d);ISchlichenm aiCT_fc_j2^ 



Borr ero et alj (120071 ) ; iBorrerol l|2007l ) ; iTritschler et all j|2007T ) ; iBellot Rubio et alj (|2004 
2003| ). These calculations demonstrate consistency between the calculated azimuthal 



variation of NCP and measurements made at low spatial resolution in visible and near 
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infrared spectral lines. Two-component inversions interpreted within the context of 
embedd ed flux tube models and applied to low spatial resolution (0.6-1 arcsec) Stokes 
data bv lBorrero et al. 1 12004 l200a . l2006h similarly were shown to be largely consistent 
with the assumed inversion (flux tube) model. 



2.2 Siphon flow and dynamic flux tube models 

The first flux tube m odel proposed to explain Evershed flows in penumbrae is the 



siphon flow model of iMever fc Schmidt! (1968). In this model, a difference in field 



strength between the two footpoints of a flux tube leads to a difference in gas pressure, 
driving a flow in the direction of the footpoint with th e hig hest field stren gth. This 
work was later followed up by IPegenhardtl dl989l. Il99lh andlThomad <ll98Sl) and th e 



model further refined in a series of p apers bv lMontesinos fc Thomasl (|l989l , ll993l , Il997f ) ; 
iThomas fc Montesinoj (1990, 1991). Given the free parameters of the calculations, the 



siphon flow model of Montesinos fc Thomas! ( 1997t ) allow consistency with the discovery 



that the Evershed flow connects to pat ches of opposite magnetic polar ity at deep layers 



near the outer boundary of a sunspot ( Westendorp Plaza et al.lll997l 'l 



Siphon flow models allow an interpretation of Evershed flows as steady flows, with 
the properties of the flow determined by assumed conditions at the footpoints of the flux 
tube and of the surrounding atmosphere. The mechanism that produces the magnetic 
field strength difference between the footpoints, needed to generate a gas p r essure 
gradient to drive the flow, is not explained by such models. Ijahn fc Schmidt] |l994) 
proposed the concept of interchange convection of magnetic flux tubes (or rather, 
sheets) as an explan a tion of the penumbral heat flux. To investigate this proposal, 
ISchlichenmaier et alJ (|l998al fbh developed a simplified ID, one-component numerical 
model of such a flux tube and studied its time evolution. In this model, a flux tube 
initially in contact with the magnetopause (the outer boundary of the sunspot) is 
heated radiatively by the external field-free gas. Its subsequent evolution is driven by 
the buoyancy of the flux tube and the superadiabatic stratification of the surrounding 
penumbra atmosphere, assumed to have properties unaffected by the flux tube. At the 
surface, radiative cooling of the tube causes it to loose buoyancy such that its upper 
part settles at a height of about 100 km above the photosphere. A gas pressure gradient, 
driving the Evershed flo w, develops along the tube f rom downstream radiative cooling. 



Later simulations bv lSchlichenmaierl (2OO2I I2003T ) with reduced numerical viscosities 



show a similar initial behavior of the flux tube. However, near the surface, the flux tube 
subsequently develops standing waves downstream from the footpoint with downflows 
diving down into the convectively unstable layers beneath the surface. The crests of this 
oscillating flux tube remain visible above the surface and show an inward migration 
in the inner part of the umbra and an outward migration in the outer penumbra 
and outside the penumbra. This behavior is similar to that of observed penumbral 
grains in the inner and outer penumbra and moving magnetic features outside the 
penumbra. The discovery of small-scale bipolar magnetic features propagating from 
the mid penumbra to outside the penumbra, where they become moving magnetic 
featur es, is consistent with the 'sea serpe nt' behavior of Schlichenmaier's moving flux 
tubes |Sainz Dalda fc Bellot Rubiolboosl ). 



A problem, investigated by several authors, e.g.. lSolanki et ahl ( 1993h ; lltezaei et al.l 



(2006), is the large radial mass flux of the Evershed flow inside the penumbra. Only 



part of this flow appears to continue in the magnetic canopy above the quiet sun 
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photosphere outside the spot. To explain this, most of the Evershed flow must submerge 
close to the outer boundary of the penumbra. The moving tube model si mulations 
show such downflows within the penumbra. How ever, Thomas] (12005 . 2006:') objected 
that the undulations seen in the 'sea serpents' of ISchlichenmaier (2002) should occur 
preferentially in the horizontal plane and hence can explain neither moving grains nor 
convective downflows. 

The question of a heating mechanism to explain the penumbral radiative heat flux 
was investigated bv lSchlichenmaier fc Solankil (|2003h . Based on estimates of the radia- 
tive cooling time and the time span of successive emergences of flux tube s, they con- 
cluded that interchange convection cannot provide the needed energy flux. IWeiss et al.l 
J2004h also argued against interchange convection on the basis that long loops of mag- 
netic field connecting to a distant active region cannot possibly interchange with hori- 
zontal fields carrying Evershed flows. 

The conclusion of lSchlichenmaier fc Solankil ( 20031 ) was that upflows along the mag- 
netic flux tubes can explain the penumbral brightness, but only if the flux tube sub- 
merges again within a distance of 1000-2000 km from their footpoint. The upflow in 
a narrow tube cannot supply the radiative energy losses over a distance correspond- 
ing to the entire radial extent of a penumbra unless it submerges and is re-heated. 
Such re-heating does not solve the energy flux problem, however, since it relies on a 
(convective) mechanism to transport the heat to the bottom of the flux tube. Never- 
theles s, the discovery of field lines r eturning to the penumbra and associated down- 
flows | W gjtgjidgrrj^laza^tal. [l997l ) was considered as support for this explanation 
( Schlichenmaier fc Solankill2003l ). ~ 



2.3 Convection and downward pumping of magnetic flux 



The siphon flow model of Montesinos and Thomas is unrelated to any convection pro- 
cess operating in the penumbra. This model represents a stationary sol ution that can - 
not explain time dependent behavior such as moving penumbral grains |Thomasll200il ). 
These grains are instead interpreted as originating from a moving convective pattern 
in the brighter parts of the penumbra (|Weissll2002l , l2006al rbl) . Whereas the moving tube 
simulations show localized downflows inside and outside the penumbra, the arched flux 
tubes of the siphon flow model re quire a mechanism to submerge and hold down their 
out er parts to su s tain e quilibrium ( Montesinos fc ThomasllT997l ) . Thomas et al. (2002a) 
and I Weiss et al. (2004) proposed that this submergence of the flux tubes occurs as the 
result of downward pumping by convection outside the sunspot. They even took this 
proposal one step further and proposed that this downward pumping is the origin of 
the filamentary structure of the penumbra. In this view, the salient features of penum- 
brae: their filamentary structures, the strong variations in magnetic field inclination 
across filaments and the Evershed flows, are to a large extent explained by what hap- 
pens outside the sunspot. Magnetic fields in brigh t and dark filament s are distinct and 
canno t be interchanged jThomas fc Weissl |2004| ). To support this, iThomas fc Weissl 
( 20041 ); IWeiss et al] (|2004l ); IWeissI |2006bl ) refer to X-ray observations and TRACE 
images showing loops extending over great distances across the Sun. We believe that 
their description and connection to the interlocking comb structure of the penumbra 
is misleading. Virtually all information about large fluctuations in the magnetic field 
inclination within the penumbra comes from spectral li nes formed withi n a few hun- 
dred km above the photosphere. The images referred to ( Sams et al.lll992l ) do not have 



5 



the spatial resolution needed to separate X-ray loops (interpreted to outline field lines) 
from bright and dark filaments. As far as the author knows, there are no observations 
that allow us to conclude that azimuthal variations in the magnetic field inclination as- 
sociated with filamentary structures correspond to field lines that are widely separated 
also far away from the penumbral photosphere. We have argued 1 Spruit fc Scharmerl 
2006) that the strong inferred variations in the magnetic field inclination within the 
first one or two hundred km above the penumbral photosphere, seen even in the in- 
ner penumbra, cannot be explained by a mechanism operating outside the sunspot. 
Instead, these strong variations suggest a local mechanism at work. The bright and 
dark filaments are not distinct and they can interchange. Furthermore, Evershed flows 
are associated with field lines that only locally and during a limited time are nearly 
horizontal. 



3 Limitations and problems of flux tube interpretations 

The success of the uncombed penumbra model i Solanki fc Montavon 19931) in explain- 



Ihe success ol the uncombed penumbra model (bolanki &: Montavon 199 6) m explain- 
ing observed Stokes spectra and ne t circular polarization (NCP) JBorrero et al.l 120071 : 
lBorrerol[2007l : iTritschler et al.ll2007l : iBellot Rubio et alJl2004 12003 ) is unquestionable. 
Furthermore, the moving tube simulations of lSchlichenmaierl (j2002j) make excellent con- 
tact with the uncombed penum bra model. It demonstrate s consistency with observed 
strong upflows in bright grains dRimmele fc Marinoll2006l) and the beha vior of bright 



grains and moving magnetic features ( Sainz Dalda fc Bellot R ubio 2008). It is hardly 



surprising that the ability of these models to explain azimuthal and line-of-sight (LOS) 
gradients of the penumbral magnetic field and Evershed flows was deemed successful. 

However, a fundamental problem remains. Whereas the embedded flux tubes sim- 
ulated by Schlichenmaier are consistent with observations i n many respects, such flux 



tubes present problems in explaining penumbral heating ( Spruit fc Scharmerl [2006: 
IScharmer fc Spruit! 120061 ). As discussed above, a horizontal flux tube is likely to heat 
the penumbra over a radial distance not much more than 1000 km ( Solanki fc Montavon! 



Il993l : ISchlichenmaie J2003I ) . which is typical of a penumbral grain rather than a penum- 
bral filament. Even over such a short distance, radiative cooling of the flow leads to 
significant temperature and brightness gradients a long the flux tube unless th ere is a 
separate source of heating below the flow channel dSchlichenmaier et al.|[l999l ). Spruit 
and Scharmer therefore argued that the presence of flux tubes covering a large fraction 
of the penumbral surface would constitute a hindrance for heating of the penumbra. 
They also pointed out that the existence of elevated flux tubes extending up to a few 
hundred km above the penumbral photosphere correspond to unlikely perturbations 
in a magnetic field so dominant already at this height that it must be expected to be 
nearly potential. Furthermore, the moving tube simulations represent a highly idealized 
model that cannot not be expected to be more than a coarse representation of reality. 
In particular: 

— The existence of the flux tube is an assumption in the model 

— The model is 1-dimensional and the flux tube assumed to be 'thin' (see below) 

— The simulations correspond to a 1-component model with the properties of the 
background atmosphere unaffected by the evolution of the flux tube 

— Only a single flux tube is simulated. The influence of neighboring flux tubes is not 
accounted for. 
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— The curvature forces of the surrounding magnetic field are ignored and its influence 
is reduced to a scalar magnetic pressure, similar to a gas pressure. 

The same objections apply to the siphon flow models discussed in previous sections. 

In view of these short-comings, it is remarkable that the moving tube simulations 
appears to capture important properties of penumbral dynamics. This is further dis- 
cussed in Sect. 4.4. 

Magnet ostatic flux tu be models including forces from a surrounding potential mag- 
netic field (jBorrerd 120071 ) demonstrate the difficulties of embedded flux tube configu- 
rations. Prescribing a specific (circular) cross section for the flux tube corresponds to 
an overconstrained problem such that not only the gas pressure but also the temper- 
ature and density within the flux tube are given by force balance alone. There is thus 
no room for an energy equation with this type of models. Furthermore, equilibrium is 
not possible with a purely potential magnetic field inside the flux tube. In the models 
shown there is an azimuthal component, corresponding to a volume current aligned 
with the flux tube, in addition to the radial field component. The bottom part of the 
flux tube is nearly evacuated whereas the top part is denser than the surroundings in 
order to balance the magnetic forces at the top and bottom, stretching and flattening 
the flux tube. These problems originate from the surrounding magnetic field wrapping 
around the flux tube and cannot be resolved by making the flux tube thinner. 

In attempt t o understand the temperature stru cture and energy balance of penum- 
bral flux tubes, iRuiz Cobo fc Bellot Rubiol {2008) developed a model for a flux tube 
with a weak magnetic field aligned with a homogeneous magnetic field along the flux 
tube axis. An objection against this model is that it suffers from a lack of consistency 
as regards force balance, which is implemented in a way that is equivalent to ignoring 
the vector properties of the surrou nding magnetic field. The origin of this problem is 
the same as that of the models of iBorrerol 1 20071 ): Prescribing the shape of the flux 
tube cross section is in general incompatible with either force balance or an energy 
equation. 



3.1 Ambiguities of interpretations based on inversions 

A major obstacle to understanding penumbra fine structure has been the lack of ade- 
quate spatial resolution in observed polarized and unpolarized spectra. In spite of suc- 
cessful adaptive optics systems operating on major solar telescopes, it has not been pos- 
sible to reach the diffraction limit with the long integration times needed for such data 
with adequate signal-to-noise. (However, by combining many short exposure frames 
and using image restoratio n techniques, filter-based systems allow near diffraction 
limited spectropolarimetry (|van Noort fc Rouppe van der Voortlliooj ; IScharmer et aJj 



2008a).) The exception is observations in th e near infrared, for exampl e those made 
with the Tenerife Infrared Polarimeter, TIP ( Martinez Pillet et al.lll999h on the Ger- 



man Vacuum Tower Telescope (VTT) at wavelengths around 1.5 /i. At that wavelength, 
the diffraction limited resolution is about 0.6 arc sec with the VTT. 

At a spatial resolution of 0.6 arc sec or worse, penumbral fine structure is not 
adequately resolved. To compensate for this and in order to test the vali dity of the em- 
bedd e d flux tube (uncombed penumbra) m odel, various observers, e.g., ( Borrero et al.l 



|2004| . 120061 ; iBellot Rubio et all 120031 . |2004| ) have implemented two-component inver 
sion techniques to interpret Stokes spectra. These investigations show that it is in- 
deed possible to fit the data with the assumed (highly idealized) representations of 
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Fig. 1 Left: Magnetostatic flux tube model of [Borrcro (2007) with dashed curves showing 
(top to bottom) optical depths r of 0.1,1 and 10. Note that only the upper part of the flux 
tube is located above the photosphere and that nearly the entire flux tube is located below 
t = 0.1. Right: Transverse field lines of a similar (but not identical) model jBorrero et al.l 
2007). Note the similarity of the magnetic field configuration above the flux tube in this model 
and that of the convective gap models, shown in Fig|4] 



flux tubes. F orward calculat i ons of Stokes spectra dMartmez Pilletl feoOO) and unpolar- 
ized spectra |Rimmelelll995l ; lBellot Rubio et al.ll2006r i based on flux tube models also 
demonstrate consistency. In some cases, it was demonstrated also that the observed 
data could equally well be reproduced wi t h flux tube representations and models with 
smooth gradients dMartmez Pillet] |2000| ; iRouppe van der Voortl 120021 ; I Borrero et al.l 
120041 : iBellot Rubio et al.ll200fj ). This ambiguity is a consequence of the width of the 
radiative transfer response function, smearing out the effects of discontinuities in the 
observed (Stokes) spectra. The interpretations of Stokes spectra clearly show com- 
patibility with flux tube interpretations. However, the simplicity of the implemented 
inversion models and the use of two components to represent observational data of 
penumbral fine structure at inadequate spatial resolution adds to these uncertainties 
to the extent that we are justified in questioning whether a description in terms of 
embedded flux tubes is an adequate representation of penumbra fine structure. 



4 Convective origin of penumbral filaments 

A n alternative explanation t o understanding penumbra fine structure was proposed 
bv lSpruit fc Scharmerl (|2006h - The filamentary structure is explained by convection in 



radially aligned (nearly) field-free gaps just below the visible surface. Such intrusions 
unavoidably lead to strong variations in the magnetic field strength and inclination 
above the gaps, but these variations are fully consistent with even a simple poten- 
tial magnetic field co nfiguration. The model explains dark cores in bright filaments 



( Scharmer et~ai1l2002l ), seen in the inner and mid penumbra (c.f., Fig.s[2]and[3]l as an 



indicator of strong field strength variations across filaments, leading to a strongly vary- 
ing Wilson depression. The Eversh ed flow is in this mode l identical to the horizontal 
flow component of the convection 1 Scharmer et alj|2008bl ). We explain this model in 
more details in the following. 



Fig. 2 S unspot located at a heliocentric distance of 20 deg, observed with the SST on 12 
Sep 2006 (va n Noort fc Rouppe van der Voortj [2008f) . The image shows the Stokes I intensity, 
averaged over the blue and red wings of the 6302 iron line. Note the dark cores, clearly visible 
in the inner penumbra. 



The large inclination gradients are a consequence of magnetic fields being diver- 
gence free (V'B = 0): field lines cannot disappear at the top of a gap (or flux tube) 
but must bend around it. To estimate the characteristic vert ical scale H for these incl i- 
nation variations, we can assume a potential magnetic field I Spruit fc Scharmerlbood ) . 
For filaments separated by a distance L, this gives H « L/2n. With a typical sepa- 
ration of 1" between filaments, this corresponds to a vertical height scale of 120 km. 
Simple magnetostatic models for such configurations, based on identical temperature 
variations with height for the two components, show distinct differences between the 
inner penumbra, where the magnetic field is more vertical and stronger than in the 
outer penumbra (jScharmer fc Spruidl2006l ). In the inner penumbra, the magnetic field 
is cusp-shaped above the gaps and associated with a large (fa 200-300 km) Wilson 
depression relative to that of the gaps. Even when the temperature is the same inside 
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Fig. 3 The same sunspot as shown in Fig[2](van Noort & Rouppc van dcr Voort 2008). The 
image shows the difference between Stokes V (circularly polarized light), recorded in the blue 
and red wings of the 6302 iron line. This serves as a proxy for the line-of-sight component 
of the magnetic field vector. Note the strongly reduced polarization signal at locations of the 
dark cores, suggesting strongly reduced field strength at these locations. 



and outside the gaps, the strong Wilson depression leads to an observed brightness 
that is lower above the gaps than between the gap s. The dark-cored filam ents discov- 
ered with the Swedish 1-m Solar Telescope (SST) l|Scharmer et all 2002). seen in the 



inner and mid penumbra, are thus explained by a combination of increased opacity 
associated with a strongly reduced fi eld strength and an overall drop of temperature 



with height (Spr uit fc Scharmerll2006D . In the outer penumbra, the Wilson depression 



is only on the order of 50 km in these models. With such a small Wilson depression, 
intensity variations from 'global' vertical temperature gradients cannot be expected to 
completely dominate over local horizontal and vertical gradients associated with details 
of the heating and cooling of convecting gas. In the outer penumbra we therefore do not 
expect the same kind of relation between filament brightness and field strength as for 
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Fig. 4 Magnetostatic convective gap models C Scharmer fc SpruitJ |2006|) for the inner (left) 
and outer (right) penumbra. The dashed curves show the height at which the gas pressure 
in the magnetic component equals the photospheric gas pressure in the field— free component. 
This serves as proxy for the r = 1 surface and leads to Wilson depressions on the order of 
200 km for the inner and 50 km for the outer penumbra. Note the similarity of the magnetic 
field configuration for the outer penumbra above the 'photosphere' in this model and the flux 
tube model, shown in Fig[T] 



the inner penumbra. The magnetostatic models are therefore in qualitative agreement 
with the absence of dark cores in the outer penumbra. 

The potential magnetic field configurations associated with two magnetostatic mod- 
els are shown in FigU The upper figure shows the calculated field lines for a weak 
(1000 G) nearly horizontal (average inclination 75 deg) magnetic field and the lower 
figure for a stronger (1800 G) and more vertical (average inclination 45 deg). Also shown 
are the shapes of the gap and the height at which the gas pressure between the gaps is 
equal to that of the field- free component at z — (horizontal dashed lines) . This serves 
as a proxy for the height at which the continuum optical depth is equal to unity. As is 
clear from the figure, the magnetic field configurations above the gaps are associated 
with strong gradients. The simple potential field model thus explains large magnetic 
field inclination variations above the penumbral photosphere without invoking forces 
in these layers. The associated current sheet is located at and below the photosphere, 
where the gas pressure is much higher than a few hundred km above the p hotosphere. 
This is in contrast to the uncombed model (|Solanki fc Montavonl [l993t ). where the 
strong gradient in the magnetic field is a direct consequence of a local perturbation in 
the form av an embedded flux tube located above the photosphere. The current sheets 
associated with such flux tubes are difficult to comb ine with magnetostatic eq uilibrium 
because of the lower gas pressure at these height s dSpruit fc Scharmeill200a ). Indeed, 
most of flux tube in the magnetostatic model of iBorrerol §200$ ) is buried below the 
photosphere and only about 130 km protrudes above the surrounding photosphere. It 
seems very difficult, if not impossible, to construct similar models for flux tubes located 
entirely above the penumbral photosphere. 

The convective gap model thus explains strong magnetic field gradients above the 
penumbral photosphere as a necessary consequence of potent ial fields. This model elim- 
inates the problems o f larg e curvature forces discussed by ISanchez Almeida fc Lite's] 
(|l992l ); ISolank7 et al.l Jl994|). constituting a corner-st one argument in favor of the un- 
combed penumbra model ( Solanki fc Montavonlfl993T ). The convective gap model also 
predicts configurations for the inner and outer penumbra that are quite different. In the 
inner penumbra, the field is cusp-shaped and associated with a large Wilson depression, 
large field strength fluctuations but relatively sm all fluctuations in inclination. In the 
outer penumbra, the magnetic field is spine- like (jLites et al.lll993T) . with a small Wil- 
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Fig. 5 Schematic drawing explaining the brightness of faculae observed near the limb 
iKeller et alj 12004ft . The low gas pressure associated with the strong magnetic field of the 
flux concentration makes it essentially transparent such that we can see the hot granular wall 
of the surrounding denser field— free gas. The optical depth unity surface is strongly tilted with 
respect to the horizontal by the shape of the flux tube magnetic field such that we see deeper in 
to the hot convecting granules close to the limb than at disk center, explaining the brightness. 
A similar situation occurs with penumbral filaments associated with strong variations in field 
strength (see FigU). Another similarity between flux tubes and faculae is that the radiative 
cooling of the surrounding field— free gas through the facula leads to a convective flow pattern 
that is downward adjacent to the facula and that this downflow is observable on the limb-side 
but hidden from view on the center-side of the facula. In a similar way, strong Wilson depres- 
sions associated with narrow penumbral filaments lead to obscuration of the limb— sides of the 
filaments already for small heliocentric distance. Convection in such filaments are also affected 
by radiative cooling through the surrounding magnetic gas such that we expect upflows in the 
middle of the filaments and downfiows at their sides. 



son depression, and small field strength fluctuations above the photosphere but with 
large inclination variations. These qualitative differe nces between the inner a nd outer 
penumbra are in good agreement with observations (jScharmer fc Spruitll2006| ). 

The overturning convective flow patterns associated with the gaps are predicted 
to be up ward in the middle and d ownflow along the boundaries to the magnetic com- 
ponents (jScharmer fc Spruitll2006l). Added to this flo w pattern is a radially outward 
(Eversh ed) flow, explained bv l Scharmer et al.l ( 2008bl ) on the basis of 3D MHD simu- 
lations ( Heinemann et al.ll2007l ) as being identical to the horizontal component of this 
convection. In our model, the dark cores of the penumbral filaments correspond to 
locations of convective upflows, in contradiction with what we expect from field-free 
convection. As explained above, the strong fluctuations in field strength across the fila- 
mentary structures in the inner penumbra lead to a correlation between brightness and 
field strength, such that we see deeper in to the hotter gas where the field strength is 
high. The mechanism for producing the dark penumbral cores i s direct l y rela ted to the 
mechanism that produces bright facul ae, proposed initially bvlSpruitl lll976h and well 
estab lished by numerical simulations (Ke ller et al.ll20o3; ICarlsson et alj|2004l ; Isteiner] 
2005) and high-resolution SST observations ( Lites et al.ll2004T ). In the case of facu- 
lae, convection occurs in the surrounding darker and field-free photosphere, whereas 
convection is inhibited or strongly suppressed in the brighter and magnetic faculae. 

The connection between dark-cored filaments and faculae can be carried further, 
c.f., Fig. [5] The brightness of faculae near the limb is explained by the r = 1 sur- 
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face being strongly inclined to the hori zontal at the limbside of faculae l|Sprurdll976l : 
iKeller et alJ 120041 : ICarlsson et"ai1l2004l ). Looking through the nearly transparent gas 
within the strong magnetic field, we see deeper into the surrounding hot convecting 
gas close to the limb than at sun center. In a similar way, we see deeper into the 
convecting parts of the penumbral filaments when viewin g sunspots well away from 
disk center and at +/- 90 degrees from the symmetry line i Scharmer fc Spruitll2006h . 
The convective gap model leads us to interpret the absence of observational evidence 
for penumbral convection as not only due to lack of spatial resolution: there is also a 
difficulty of seeing deep enough into the filaments to be able to observe the convection 
for sunspots close to disk center. In addition, there is the confusion from the brightness 
- field strength correlation already discussed above. 



4.1 Limitations of the convective gap model 



The magnetostatic gap models discussed above ( Scharmer fc Spruidl2006l ) predict a 



grad ual transition fro m cusp-shaped magnetic fields in the inner penumbra to spine- 
like i Lites et al.lll993T) magnetic field configurations in the outer penumbra, in good 



agreement with observations. However, these simple models assume a perfectly field- 
free gap and do not include an energy equation, nor are the forces associated with 
the convective flows included. Details of flows and their interactions with the magnetic 
field cannot be explained with this simple model. The explanation of dark cores re- 
lies on opacity effects that are obvious only for the model corresponding to the inner 
penumbra. We expect these results to be relatively robust. In the outer penumbra, pre- 
dictions about filamentary brightness are more difficult without an energy equation. 
Furthermore, observations show strong Evershed flows in the outer penumbra where 
the magnetic field is weaker. This combines to making the kinetic energy density pv 2 /2 
of comparable magnitude to the magnetic energy density B 2 /2po such that we expect 
relatively strong effects from the flow on the magnetic field. Whereas the magneto- 
static gap models show good overall agreement with observed properties of penum- 
bral magnetic fields, we cannot expect detailed agreement between the gap model and 
observations also in the outer penumbra. Only more accurate models and numerical 
simulations can provide this. 



4.2 Support from observations 



The interpretation of lig ht bridges as essentially field-fr ee gaps dividing the umbra of 
a sunspot in two parts |Lekalll997l ; I Jurcak et ahl l2006h does not seem controversial. 
3D MHD simulations of such structures ( Nordlundl 20061 : iHeinemanrj 20061) reproduce 
observed dark lanes running along the center of such structures ( Lites et ahl l2004h 
and demonstrate that the origin of this dark structure is the same as proposed for 
the convective gap model. Of considerable importance is therefore that dark-cored 
light bridg e structures occ a sionally show smooth transitions to dark-cored penumbral 
filaments ( Langhansl 120061 : IScharmer et al.1 120071 ), strongly suggesting simil ar origin. 
Upflows in light-bridge dark lanes and the dark cores of penumbral filaments ( Rimmelg 
2008) suggest a common interpretation in terms of convection, but with evidence for 
horizontal flows at greater heights also fitting a flux tube interpretation. Connections 
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Fig. 6 Examp les of space-time sli ces across filaments in the inner penumbra observed with 
SOT/Hinodc (Ichimoto et al. 2007). The disk center direction is indicated with an arrow and 
labeled with the corresponding heliocentric distance. The upper image of each pair shows the 
surroundings of the cut across the filaments (dashed line) analysed in the space-time slices 
(lower image). 



of dark-cored penumbral filaments to periphera l umbra l dots and dark cores in light 
bridges have been reported also bv lBharti et al. | |2007bl ). 

Several recent papers report evidence for c onvection in umbral dots dBharti et al.l 
l2007al ; iRimmelel I200I iRiethmiiller" et al.l l2008). Although this provides no direct ev- 
idence for penumbral convection, the direct connection of peripheral umbr al dots to 
dark-cored filaments ( Langhansl I2OO6 ; Lan ghans et alj|2007l ; lRimmele| [2008) provides 
'circumstantial' evidence for this interpretation. We ca ution, however, that the um bral 
dots observed are much larger than those simulated bv lSchiissler fc Voglerl |2006l ) and 
in some cases resemble granular intrusions. 

In a paper of fundamental importance, llchimoto et alj (|2007m find strong evidence 
for overturning penumbral convection by analyzing continuum images and spectra ob- 
served with SOT on Hinode. Space-time plots created along lines crossing filaments in 
the inner penumbra, at directions +/- 90 deg from the line connecting the center of 
the solar disk and the center of the sunspot, displayed twisted ropelike structures. The 
'twisting' motion was consistently in the direction toward sun center for both sides of 
the spot and irrespective of whether the spot was East or West of the meridian, (c.f, 
Fig©- The apparent twist observed is therefore neither an actual twist nor a helical 
motion of individual filaments but must be a viewing angle effect. The interpretation 
(|lchimoto et al.l l2007t) is of upflows of overturning convection, viewed from the side. 
With the limb-side part of these filaments hidden from view, such flows will always 
appear to be in the direction of sun center direction for spots observed away from 
disk center. This is in perfect a greement with the predict ions of the magnetostatic gap 
model for the inner penumbra (jScharmer fc SpruitJbOQfJ) . Here the Wilson depression 
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between the gaps is predicted to be so large that the limbside part of the filaments is 
invisible already for disk center distances in the range 21-35 deg. The spot observed 
by Ichimoto et al. was located at a disk cente r distance of 31-48 deg. Th e flow pattern 
observed is consistent with that predic ted by Scharmer fc Spruit] ( 2006h . 

Based on high-resolution SST data. Izakharov et al. ( 2008h recently inferred similar 
evidence of convective flows for a sunspot loc ated 40 deg from disk center. Surprisingly, 
this is interpreted in terms of convective rolls (|D anielsonll 1 96 lh . However, the horizontal 
flow toward the center of the filament at the bottom of such a roll, needed to verify its 
existence, is unobservable. Furthermore, the life times of penumbral filaments are on 
the order of 1 hour or more l Langhans et al.ll2007l ). To sustain the radiative output over 
such a long time, the convective upflow must persist to depths much larger than a few 
hundred km. Finally, the observations of llchimoto et al. I (|2007h were made in the inner 
penumbra, where the magnetic field has a strong vertical component, whereas roll-like 
convection is expected to be primarily associated with more horizontal magnetic field. 



4.3 Support from 3D MHD simulations 



In contrast to what is to what is the case for faculae ( Keller et al.l l2004; Car lsson et alj 



20041 ) and umbral dots (Schiis sler fc Vogle rl2006h . realistic simulations of entire sunspots 



have not yet been feasible. This is partly due to the difficulties of thermally relaxing 
such a deep structure and maintaining its stability but mostly due to the huge range 
of scales associated with a fully developed sunspot. The first attempts to carry out 3D 
MHD simulations with radiative energy tra nsfer of fine structure in a small sunspot 
were ca rried out by iHeinemann et aD (]2007h and further discussed by IScharmer et al.l 
(|2008b| y Recently, simulations using a similar setu p and grid separation , but with a 
ich larger computational box, were carried out bv lRempel et al | |2008h . A synthetic 



muc 

continuum image calculated from these simulations is shown in Fig. [7] The approach 
taken in both simulations is to reduce the computational effort by using a rectangular 
computational box containing only a small 'azimuthal' slice of a sunspot. The overall 
results of these simulations obtained with two independent codes are quite similar, 
although differing strongly in the length of the penumbral filaments: 

— The origin of filamentary structures is overturning convection and the dark cores 
are caused by a locally ele vated r = 1 surface, supporting the convective gap model 
( Spruit fc Scharmer1l2006l l. 



— The convection occurs in deep gaps, up to about 2 Mm (Re mpel et al.|[2008l ) with 
strongly reduced field strength. 

— The simulations show horizontal outflows, similar to Evershed flows but with smaller 
velocities, peaking near optical depth unity and associated with locally strongly in- 
clined fields. 

— The bright heads of the penumbra filaments show inward propagation and strong 
upflows 

— The simulations show moving magnetic features (MMF's) and moat flow in the 
surrounding photosphere 

For overturning convection to be efficient, the gas needs to stay near the surface 
for a significant amount of time in order to give it time to cool radiatively. At the 
same time, it needs to move away from its upflow point in order to allow more gas to 
flow up. Horizontal flows are thus essential components of overturning convection. The 
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Fig. 7 Synthetic continuum ima ge calculated from 3 D MHD simulations, showing a sunspot 
and its surrounding photosphere (Rempel ct al. 2008). In the umbral part of the spot is seen 
umbral dots with dark cores, in the penumbral part filaments that reach lengths of up to 
2-3 Mm. 



Evcrshed flow is identi fied as being identical t o the horizontal flow component of this 
penumbral convection ( Scharmer et alJl2008bl ). 

The simulations reproduce fundamental properties of observed penumbrae. This 
gives us confidence in concluding that the simulations constitute good representations 
of penumbral dynamics and energy balance. However, neither of these penumbra sim- 
ulations show nearly horizontal magnetic fields in the outer penumbra, as observed, so 
something is missing in the simulations. Other features that do not appear realistic are 
elongated structures in the outer penumbra that appear intermediate to granules and 
filaments and a less distinct boundary to the umbra than observed. A complex question 
concerns the field strengths in the gaps. In the simulation s these are on the order of 
700-1000 G. In the simulations of iHeinemann et al.l (|2007l ). this was interpreted as a 
consequence of numerical diffusivities, allowing overturning convective flows to cross 
field lines. Similar processes take place in the simulations of Rempel et al. (Rempel, 
private communication). The existence of strong magnetic fields associated with over- 
turning convection in penumbrae thus rely on turbul ent magnetic diffusio n at small 
scales. I n contrast, umbral dots in the simulations of iRempel et al.l (|2008i ) as well as 
those of lSchiissler fc Voglerl 1 20061 ) appear as nearly field-free plumes near the surface. 

The simulations of Rempel et al. show som e indication of flow patterns that are 
reminiscent of roll convection ijDanielsonl 1 1 96 1 ) . but this interpretation is uncertain. 
The convective energy transport is mainly provided by deep upflow and downflow 
plumes. An intriguing result of both simulations is that the convective gaps do not 
extend to the field- free atmosphere below the sunspot, suggesting that sunspots radiate 
energy that is contained within the magnetic field initially. A difference between the 
simulations of Rempel et al. and Heinemann et al. is that the former simulations show 
radial inflows adjacent to the outflows near the surface whereas the latter simulations 
show evidence of such return flows only deeper below the surface. 

A perhaps relevant result for understanding penumbral magnetic fields comes from 
simulations of quiet sun emerging flux. In these simulations, the expulsion of magnetic 
flux takes place in the horizontal dire ction by horizontal flows but also in the vertical 
direction by overshooting convection (Stci ner et al.| [2008). The weak horizontal fields 
of such emerging flux are expelled to heights of about 500 km, where they are difficult 
to observe in photospheric lines. Such expelled magnetic fields are limited to a height 
at which the gas pressure below the expelled magnetic field is roughly equal to B 2 /2[iq. 
If a similar expulsion mechanism operates in the penumbra, where the field strength is 
much higher, any expelled horizontal field must be located much closer to the penum- 
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bral photosphere than for the emerging flux discussed bv ISteiner et al.l 1 2008h . Such 
magnetized gas would then show strong polarization signatures reminiscent of flux 
tubes. 



4.4 Connections to flux tube models 



In spite of its ID rep resentation and the failure to c onfirm interchange co nvection, the 
moving tube model ( Schlichenmaier et al.l Il998al lbl; ISchlichenmaierl 120021 ) has connec- 



tions to convective processes. The upflow within the flux tube is driven by the superadi- 
abatic stratification of the external atmosphere, similar to that of a field- free convective 
upflow. The flux tube evolution is similarly driven by the superadiabatic stratification 
of the external atmosphere, even though the complete cycle of inward/outward move- 
ment and heating/cooling of the flux tubes does not take place in the simulations: a flux 
tube initially located along the magnetopause moves toward the umbra and remains 
there. These ID s imulations cannot sho w convective flow patterns. However, the 'sea 
serpent' behavior (|Schlichenmaierll2002f) . with upflows and downflows along the length 
of the tube, can be in terpreted as a ID representation of 3D con vection. In the 3D 
penumbra simulations ( Heinemann et al1l2007l : lllempel et alJl20oil ). the upflows occur 



at the centers of the gaps and the downflows on either side of the upflows, displaced 
both the in the azimuthal and radially outward directions. In the moving tube model, 
such upflows and downflows can only be spatially separated in the radial direction. 
Both models have in common a convective upflow and radiative cooling driving an 
outflow away from the center of the spot. In both moving tube and 3D penumbra 
simulations, the outflow peaks in a thin layer near r = 1, which is where the gas cools 
most efficiently ( Scharmer et al.ll2008bl ). 



In the convective gap model, the strong magnetic field gradients above the gap are 
explained as a perturbation of a nearly potential magnetic field above the penumbra, 



intro duced by the nearly field-free gap ( Spruit fc Scharmer! 120061 ; IScharmer fc Spruitl 



This leads to a magnetic field that is cusp-shaped in the inner penumbra and 
locally nearly horizontal in the outer penumbra. Adding a weak horizontal magnetic 
inside the gap will not change this configuration significantly. Such a configuration is 
in its upper parts quite similar to a flux tube but deeper down, these two types of 
structures are very different. 

Support for the existence of embedded flux tubes based on a magnetostatic mode l 
is claimed from calculations of net circular polarization (NCP), ( Borrero et al.ll2007l ). 



However, the r = I surface of these models intersect the symmetry axis of the flux 
tube well above the center of the flux tube and the r = 0.1 surface (typical of the line 
formation height) cuts through the top of the flux tube (c.f., Fig.[T}. These calculations 
clearly are sensitive only to the upper part of the flux tube, where its magnetic field is 
similar to that of the convective gap model. 

A similar ambiguity concerns t he origin of the dark cores of penumbral filaments, ex - 
plained by convective gap models jSpruit fc Scharmerlliooj ; IScharmer fc Spruitlbood ). 
Also flux tube models with weaker field in the flux tub e than in the surroundings 
rod uce opacity effects resulting in dark-cored structures ( Ruiz Cobo fc Bellot Rubiol 
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2008). This is interpreted as support for flux tube models by the authors. However, also 
for this model, only the top of the flux tube is visible above r = 1, so this configuration 
is similar to a convective gap model in its observable parts. 
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Interpretations of highly resolved Stokes spectra SOT/Hinode show wrapping around 
structures that can similarly be interpreted either as flux tubes or as convecting gaps 
( Borrero et al J2008h . These and other recent high-resolution Stokes data do not provide 
evidence for flow c hannels and flux tubes elevat e d above the phot o sphere, as discussed i n 
some papers, e.g.. ISolanki fc Montavonl [|l993h ; iRimmelel ifliil ) : iBorrero et alJ l|2006l ). 
Based on observations of penumbra magnetic fields, interpretations in terms of con- 
vective gaps or flux tubes partly buried below the r = 1 surface thus are inherently 
ambiguous. This also serves as a reminder of the difficulties of interpreting (inade- 
quately resolved) observations in terms of unknown underlying physics. 

While this ambiguity, in our opinion, undermines arguments for the very existence 
of embedded flux tubes, it primarily suggests that magnetic field measurements are not 
likely to show a clear distinction between flux tube and convective gap models, at least 
in the outer parts of the penumbra. The distinction between the models may need to be 
based primarily on measurements of the velocity field, which is horizontal and along a 
flux tube but with added vertical and azimuthal flow components in the convective gap 
models. It is this diagnostics that so far provides the stro ngest observational evidence 
in favor of the convecting gap model ( Ichimoto et al.ll2007l ). The observational evid ence 
accumulated so far is however too scarce to be conclusive. IZakharov et al. also 
reported evidence for convective upflows in penumbral filaments, but clearly observa- 
tions of the entire sequence of upflows, horizontal flows and downflows are needed to 
fully reveal the nature of convection in penumbrae. 

A highly controversial issue is wh ether the penumbra l conve ction is essentially field- 
free or associated with kG strengths. iBellot Rubio et al.l ( 200?)), based on SOT/Hinode 
data and Milne-Eddington (ME) inversi ons, found only small var iations in field strength 
across dark-cored penumbral filaments. IScharmer et all (|2008al ). based on SST obser- 
vations and ME inversions, demonstrated that an improvement of the spatial resolu- 
tion from 0.3" to 0.15" increases measured field streng th variations over dark cores 
by approximately a factor of two. IZakharov et al.1 ( 20081 ) . based on SST data and ME 
inversions, found loc ally weaker fie l ds by approximately a factor two, associated with 
convecting filaments. lJurcak et al. I (|2007T ) using SOT/Hinode data and inversions al- 
lowing for gradients and a Gaussian perturbation in the magnetic field to represent 
flux tubes, concluded that the field strength is reduced by only 600 G at the centers of 
bright filaments (where dark cores should be located) in the inner penumbra and fur- 
thermore that this reduction in field strength occurs only close to the photosphere and 
disappears already at logr = —0.5. Based on observations and inversions, there is so 
far no support for the assumption that these gaps are nearly field-free. This conclusion 
refers to observational data from the layers above the photosphere whereas the nature 
of penumbral convection in deeper layers can ultimately only be determined from 3D 
MHD simulations. 



5 Conclusions 

We believe that there is now strong evidence to support the conclusion that penumbral 
fine struc ture should be inte r preted as the result of overturning convection, as pro- 
posed by ISpruit fc Scharmerl l|2006h . Evidence for this conclusion comes from recent 
SOT/Hinode and SST observations (jlchimoto et al.ll2007l ; [Zakharov et alj|2008j ). show- 
ing vertical flows of the right magnitude to explain the penumbral radiative heat flux . 
Recent numerical 3D MHD simulations ( Heinemann et aTll2007l ; iRempel et al.l [2008) 
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reproduce fundamental properties of observed penumbrae and confirm the convective 
origin of penumbral filaments. The simulations show that the nature of this convec- 
tion takes place in gaps with up to 2 Mm depth and tha t any roll-like convection 
( Danielsonlfl96ll;lThomas fc Weissll2004l ; iRempel et al.ll20ol ). if present, is of small im- 
portance ( Rempel et al, 20081 ). The Evershed flow is interpreted to b e identical to the 
horizontal flow component of this convection ( Scharmer et al.ll2008bh . Such horizontal 
flows are necessary in order to cool hot upflows by radiation. 

Neither observations nor simu lations lead to the conclu sion that this convection 
is nearly field-free, as suggested ( Spruit fc Scharmerl bOPril ). However, inferred field 
strengths from spectropolarimetric data are obviously limited to layers above the pho- 
tosphere, whereas simulations rely on numerical diffusivities to prevent instabilities at 
scales corresponding to the grid separation. Other uncertainties relate to the outer 
parts of penumbrae where observati ons show nearly horizont al field and even field lines 
dipping down into the photosphere IjBorrero fc Solankill2008l ) . Simulations do not show 
structures of this type. It appears likely that downward pumping of magnetic field by 
convection outside the sunspot plays a role in the outer penumbra, as proposed earlier 
( Thomas et"ai1l2002bh . but we disagree strongly with the conclusion that this explains 
the origin of the filamentary structure of the penumbra. Downward pumping by the 
convection in side the penumbra also must t ake place and this probably explains why 
observations (|Westendorp Plaza et alj 19971) show evid ence of return flux well inside 
the outer penumbral boundary ( Scharmer et al.l 1200801 ). 

Penumbral filaments have been successfully interpreted in terms of embedded flux 
tubes during a period of 15 years. While we conclude that this interpretation is mis- 
leading in terms of underlying physics, there are several reasons why this model has 
been so successful. We have shown that the opening of radially aligned gaps with nearly 
field-free convecting gas leads to a magneti c field that is much more horizontal over 
the gaps, giving the illusion of a flux tube ( Spruit fc Scharmerl bOQlj ). There are also 
other arguments for expecting nearly horizontal fields in the penumbral atmosphere: 
Horizontal cooling flows are most efficient near optical depth unity and if this gas is 
magnetized, it will aid in producing nearly horizontal magnetic fields. Also, emerging 
flux simulations relevant to the quiet sun suggests that convection can lead to flux 
expulsion in the vertical direction in addition to the horizontal direction and that this 
explains quiet sun horizontal magnetic fields above the photosphere. Clear evidence of 
such 'vertical ' flux expulsion is ho wever not seen in the penumbral part of the sunspot 
simulated bv lRempel et al.l (|2008l ). 

We emphasize the connections of the moving tube model to convective processes 
and to the radiative cooling of such flows near the photosphere in both types of models. 
The similarity of the magnetic field above convecting gaps and flux tubes add to the 
difficulties of correctly interpreting observations and distinguishing between models. 

In the embedded flux tube models, the Evershed flow is at center stage and the 
mechanism for heating the penumbra remains obscure. The new view of penumbral fine 
structure as caused by overturning convection implies that the main driver of penumbra 
fine structure is the energy flux b elow the surface and th at the Evershed flow is 'only' 
a consequence of this convection (|Scharmer et al.ll2008bl). 

We believe that siphon flow models 1 Montesinos fc Thomas! Fl997l ) are of little rel- 
evance for understanding penumbrae. These are linked to the idea that there are two 
distinct families of field lines: those associated with dark filaments and Evershed flows 
and those associa t ed with bright filamen ts connecting to distant magnetic regions 
( Weiss et alj |2004 iThomas fc Weissll2004l ). As far as we know, there is no observa- 
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tional support for this 'static' picture of penumbral magnetic fields. Observations sug- 
gest life times for penumbral filaments on t he order of one hour asso ciated with flow 
channels opening and closing continuously jRimmele fc Marinoll200d ). 3D MHD Sim- 
ulations ( Rempel et al.ll2008l ; llfeinemann et all 2007 ) confirm the transient nature of 
azi muthal variations in field s t rength and inclination. T heoretical arguments and mod- 
els (|Spruit fc Scharmeril2006l : IScharmer fc Spruitl l200d ) as well as simulations clearly 
lead to the conclusion that the large variations in inclination across filaments are local 
perturbations, caused by penumbral convection and vanishing a few hundred km above 
the penumbral photosphere. 



As regards further progress in this rapidly evolving field, we expect that even more 
realistic 3D MHD simulations in the near future will further improve our understanding 
of penumbrae, in particular as regards their outermost parts. Observed unpolarized and 
polarized spectra at the highest possible spatial resolution are needed. Of particular 
importance is such spectra giving information about the layer immediately above the 
photosphere. Emphasis should be given to analyzing data at +/- 90 de g from the sym- 
metry axis of sunspots located away from disk center, as was done by Ichimoto et alj 
( 20071 ). This is in part to allow analysis of flows perpendicular to the radial direction of 
the filaments, but also in order to see as deep into these structures as possible. Analysis 
of such data need to account for the pronounced 3D nature of these filaments, caused 
by strong azimuthal variations in the Wilson depression, as well as strong LOS varia- 
tions in the magnetic field and flow velocity. A dilemma here is that the use of inversion 
techniques with many nodes along the LOS raises questions of uniqueness and diffi- 
culties in comparing the results of such inversions with simulations. Existing 3D MHD 
simulation data allow inversion techniques to be tested with synthetic Stokes spectra 
from penumbral atmospheres, as done already with simulations of small-scale flux con - 
centrations outside sunspots (|Khomenko fc ColladoslfeoOTl ; lOrozco Suarez et al.ll2007T ). 
The effect of assuming e.g. hydrostatic equilibrium can be evaluated quantitatively. 



Presently used inversion techniques process polarized spectra pixel by pixel without 
constraining, for example, the magnetic field to be divergence- free. The requirement 
of divergence-free magnetic fields is crucial in forcing field lines to bend around con- 
verting gaps (and flux tubes) , leading to strong gradients in field strength and inclina- 
tion. With spectropolar imetric observations approaching a spatial resolution of 100 km 
( Scharmer et al.ll200 8a). which is similar to the equivalent LOS resolution achieved 
with inversion techniques using a small number of nodes, it is reasonable to enforce 
magnetic fields constrained by div(B)=0. Stray-light corrections are with most inver- 
sion techniques implemented in an ad-hoc manner pixel by pixel whereas a physical 
stray-light implementation would employ a point spread function that does not vary, 
or varies slowly, across the FOV. Micro- and macro-turbulence parameters are used as 
fudge parameters to compensate spatial smearing of unresolved structures. With im- 
proved spatial resolution, modeled LOS velocity gradients should eliminate the need for 
such parameters in the inversions. We expect future inversion techniques to develop 
as 'global' techniques, in the sense of fitting model parameters for a large number 
of connected pixels simultaneously. This will allow constraints, such as div(B)=0, and 
physical straylight models to be incorporated in a consistent manner to further enhance 
the usefulness of inversion techniques for inferring the physical state of the atmospheres 
above sunspots and other magnetic structures. 
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